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Abstract 

This report deals with the basic concepts on deducing transit times for quantum scattering: 
the stationary phase method and its relation with delay times for relativistic and non-relativistic 
tunneling particles. After reexamining the above-barrier diffusion problem, wc notice that the 
applicability of this method is constrained by several subtleties in deriving the phase time that 
describes the localization of scattered wave packets. Using a recently developed procedure - multi- 
ple wave packet decomposition - for some specifical colliding configurations, we demonstrate that 
the analytical difficulties arising when the stationary phase method is applied for obtaining phase 
(traversal) times arc all overcome. In this case, we also investigate the general relation between 
phase times and dwell times for quantum tunneling/scattering. Considering a symmetrical col- 
lision of two identical wave packets with an one-dimensional barrier, we demonstrate that these 
two distinct transit time definitions are explicitly connected. The traversal times are obtained for 
a symmetrized (two identical bosons) and an antisymmetrized (two identical fermions) quantum 
colliding configuration. Multiple wave packet decomposition shows us that the phase time (group 
delay) describes the exact position of the scattered particles and, in addition to the exact relation 
with the dwell time, leads to correct conceptual understanding of both transit time definitions. 
At last, we extend the non-relativistic formalism to the solutions for the tunneling zone of a one- 
dimcnsional electrostatic potential in the relativistic (Dirac to Klein-Gordon) wave equation where 
the incoming wave packet exhibits the possibility of being almost totally transmitted through the 
potential barrier. The conditions for the occurrence of accelerated and, eventually, superluminal 
tunneling transmission probabilities are all quantified and the problematic superluminal inter- 
pretation based on the non-relativistic tunneling dynamics is revisited. Lessons concerning the 
dynamics of relativistic tunneling and the mathematical structure of its solutions suggest revealing 
insights into mathematically analogous condensed-matter experiments using electrostatic barriers 
in single- and bi-layer graphene, for which the accelerated tunneling effect deserves a more careful 
investigation. 
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I. INTRODUCTION 

The analytical methods utilized for reproducing the wave packet collision with a po- 
tential barrier have been widely discussed in the context of (one-dimensional) scattering and 
tunneling phenomena [1-3]. Within this extensively explored scenario, the stationary phase 
method (SPM) is the simplest and the most common approximation adopted for obtain- 
ing the group velocity of a wave packet in a quantum scattering process which represents 
the collision of a particle with a square (rectangular) barrier of potential. First introduced 
to physics by Stokes and Kelvin [4], the SPM essentially enables us to parameterize some 
subtleties of several quantum phenomena, such as tunneling [1, 5, 6], resonances [7-9], 
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incidence-reflection transmission interferences [2] as well as the Hartman effect [10] and its 
superluminal traversal time interpretation [3, 11, 12]. It is also often quoted in the context 
of testing different theories for temporal quantities such as arrival, dwell and delay times 
[1, 3, 13] and the asymptotic behavior at long times [11, 14]. The SPM is the simplest and 
most usual approximation method for describing the group velocity of a wave packet in a 
quantum scattering process represented by the collision of a particle with a potential barrier 
[3, 10, 12, 15, 16, 19]. 

To implement the accurate analysis with the SPM, it is necessary to work with the 
hypothesis of adding up the Fourier modes from the stationary solution to observe the evo- 
lution of a linear wave packet on collision with an obstacle. Those interested in linear theory 
often look at stationary (frequency-domain) solutions of a single frequency, whereas those 
interested in nonlinear theory usually look at time-dependent pulses. The time-dependent 
scattering of a wave-packet by a potential barrier in the linear Schrodinger equation is in- 
teresting because it reconciles two types of wave theories. Their limitations, however, does 
not concern the domain of the method, i. e. its mathematical foundations and applicability, 
they focus on the indiscriminate use of the method for time-dependent scattering problems, 
in particular, for deriving tunneling phase times. 

The issue of tunneling time is replete with controversy and paradoxes, preponderantly, 
because many tunneling time definitions seem to predict superluminal tunneling velocities. 
Tunneling occurs when a wave impinges on a thin barrier of opaque material and some 
small amount of the wave leaks through to the other side. In particular, the Hartman 
effect [10] predicts that the tunneling time becomes independent of barrier length for thick 
enough barriers, ultimately resulting in superluminal tunneling velocities. In this scenario, 
the correct quantification of the analytical incongruities which restrict the applicability of 
the principle of stationary phase is extended to the study of tunneling where, eventually, 
the superluminal interpretation of transition times can be ruined. 

Trying to understand the Nature of the superluminal barrier tunneling has brought up a 
fruitful discussion in the literature [12, 13, 20, 21] since pulses of light and microwaves appear 
to tunnel through a barrier at speeds faster than a reference pulse moves through a vacuum 
[5, 22-25]. The experiments that promoted the controversial discussions were performed 
with a lattice of layers of transparent and opaque materials arranged so that waves of 
some frequencies are reflected (through destructive interference) but other frequencies pass 
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through the lattices in a kind of filter eSect correlated with the Hartman effect [10]. 

To find a definitive answer for the time spent by particle to penetrate a classically forbid- 
den region delimited by a potential barrier, people have introduced quantities that have the 
dimension of time and can somehow be associated with the passage of the particle through 
the barrier, i. e. the tunneling times. These proposals have led to the introduction of 
several theoretical definitions [6, 11, 12, 15, 26-33], some of which can be organized into 
three groups. (1) The first group comprises a time-dependent description in terms of wave 
packets where some features of an incident packet and the comparable features of the trans- 
mitted packet are utilized to describe a quantifiable delay as a tunneling time [1]. (2) In the 
second group the tunneling times are computed based on averages over a set of kinematical 
paths, whose distribution is supposed to describe the particle motion inside a barrier. In this 
case, Feynman paths are used like real paths to calculate an average tunneling time with 
the weighting function exp [i S x{t)/h], where S is the action associated with the path x{t) 
(where x{t) represents the Feynman paths initiated from a point on the left of the barrier 
and ending at another point on the right of it [34] ) . The Wigner distribution paths [6] , and 
the Bohm approach [35, 36] are included in this group. (3) In the third group we notice the 
introduction of a new degree of freedom, constituting a physical clock for the measurements 
of tunneling times. This group comprises the methods with a Larmor clock [26] or an os- 
cillating barrier [37] . Separately, standing on itself is the dwell time defined by the interval 
during which the incident flux has to exist and act, to provide the expected accumulated 
particle storage, inside the barrier [3]. 

There is no general agreement [12, 47] among the above deflnitions about the mean- 
ing of tunneling times (some of the proposed tunneling times are actually traversal times, 
while others seem to represent in reality only the spread of their distributions) and about 
which, if any, of them is the proper tunneling time [12] . In the context on which we intend 
to elaborate, the tunneling mechanism is embedded by theoretical constructions involving 
analytically-continuous Gaussian, or inflnite-bandwidth step pulses to examine the tunnel- 
ing process. Nevertheless, such holomorphic functions do not have a well-deflned front in a 
manner that the interpretation of the wave packet speed of propagation becomes ambiguous. 
Moreover, infinite-bandwidth signals cannot propagate through any real physical medium 
(whose transfer function is therefore finite) without pulse distortion, which also leads to 
ambiguities in determining the propagation velocity during the tunneling process. 
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In what concerns such ambiguities, some of the barrier traversal time definitions lead, 
under tunneling time conditions, to very short times, which can even become negative. 
It can precipitately induces an interpretation of violation of simple concepts of causality 
Otherwise, negative speeds do not seem to create problems with causality, since they were 
predicted both within special relativity and within quantum mechanics [32]. A possible 
explanation of the time advancements related to the negative speeds can come, in any 
case, from consideration of the very rapid spreading of the initial and transmitted wave 
packets for large momentum distribution widths. Due to the similarities between tunneling 
(quantum) packets and evanescent (classical) waves, exactly the same phenomena are to 
be expected in the case of classical barriers [77]. The existence of such negative times is 
predicted by relativity itself based on its ordinary postulates [12], and they appear to have 
been experimentally detected in many works [38, 39]. 

In all cases described by the non-relativistic (Schrodinger) dynamics [12] , the pulse (wave 
packet) that emerges from the tunneling process is greatly attenuated and front-loaded due 
to the filter effect (only the leading edge of the incident wave packet survives the tunneling 
process without being severally attenuated to the point that it cannot be detected). If one 
measures the speed by the peak of the pulse, it looks faster than the incident wave packet. 
Since the transmission probability depends analytically on the momentum component k 
(T = T{k)) the initial (incident wave) momentum distribution can be completely distorted 
by the presence of the barrier of potential. As there is no sharp beginning to a pulse, we 
cannot declare the instant of its arrival at a certain point. Thus the computation of the 
tunneling time becomes fundamentally meaningless. We could only watch the rising edge of 
the pulse and try to recognize what is arriving. 

Meanwhile, in some limiting cases of the relativistic (Klein) tunneling phenomena where 
the relativistic kinetic energy is approximately equal to the potential energy of the bar- 
rier, and mcL/h « 1, particles with mass m can pass through a potential barrier of 
width L with transmission probability T approximately equal to one. Differently from the 
non-relativistic tunneling analysis, the original momentum is kept undistorted and there is 
no filter effect. Our final aim is thus to demonstrate that the tunneling time is then com- 
puted for a completely undistorted transmitted wave packet, which legitimizes any eventual 
accelerated transmission. 

In this scenario, the first part of our study is to characterize some limitations on the 
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use of the SPM for obtaining the right position of a propagating wave packet subject to 
the most common one-dimensional potential configuration. In section II, the analysis of the 
applicabihty of the SPM is concentrated on the investigation of the non-relativistic diffusion 
of an incoming single wave packet with energy spectrum totally localized above the potential 
barrier. We are particularly interested in obtaining an analytical description of the above- 
barrier diffusion problem by pointing out some interpretive problems which concern with 
non-conservation of probability during the collision process. By recurring to a wave packet 
multiple peak decomposition [16], we can decompose the original colliding wave packet 
into multiple reflected and transmitted components by means of which the conservation of 
probabilities is recovered during the scattering process. The analytic expressions for all the 
reflected, transmitted and intermediary components are obtained and the validity of them 
is discussed in terms of an illustrative example where the analytic approximation and the 
reflects exactly the numerical results. 

In section III, we present the theoretical constructions of the physical framework where 
the study of a tunneling process is embedded. It involves analytically-continuous Gaussian 
pulses which due to its not well-defined front lead to ambiguities in the interpretation of 
the wave packet speed of propagation. Beside it, infinite bandwidth signals cannot propa- 
gate through any real physical medium (whose transfer function is therefore finite) without 
pulse distortion, which also leads to ambiguities in determining propagation velocity during 
the tunneling process. Moreover, some of the barrier transposing time definitions lead, in 
tunneling time conditions, to very short times which can even become negative where they 
may seem to contradict simple concepts of causality. 

To partially overcome these incompatibilities, in section IV, we discuss a theoretical 
exercise involving a symmetrical collision with a one-dimensional square potential where 
we recompose the scattered wave packet by summing the amplitudes of the refiected and 
transmitted waves in the scope of what we classify as a multiple peak decomposition analysis 
[16], where the conditions for the SPM apphcability are totally recovered. 

The general and explicit relation between phase times and dwell times for quantum 
tunneling or scattering is investigated in section V. For the quoted symmetrical colliding 
configuration, the phase time gives the exact position of the center of mass of each symmet- 
rically scattered wave packet. Consequently, we can have a realistic idea of the magnitude 
the dwell time for the case in which the energy of the incident particle is smaller than the 
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barrier potential energy (tunneling configuration). In spite of a theoretical focus, the results 
here obtained apply only to such configurations which should deserve further attention by ex- 
perimenter, while the existing experiments report the inaccurate results on non-symmetrical 

configurations. 

In section VI we compute the tunneling transmission probabilities and the respective 
delay-times for wave packets when a relativistic dynamic is taken into account [17, 18]. The 
transmission probabilities, the phase times and the dwell times for the relativistic dynam- 
ics are all obtained and the conditions for the occurrence of accelerated and, eventually, 
superluminal tunneling transmission probabilities are specified. 

We draw our conclusions in section VII pointing out suggestive insights into condensed- 
matter experiments using electrostatic barriers in single- and bi-layer graphene, for which 
the accelerated tunneling effect deserves a more careful investigation. 

II. THE STATIONARY PHASE METHOD IN THE SCATTERING PROBLEM 

The principle of stationary phase to describe the group velocity is well-estabhshed, and 
was employed by Sommerfeld and Brillouin in their early studies of wave propagation and 
group velocity [40]. The method has in time become a standard tool for several theoretical 
applications [41] adopted not only by physicists but also by biologists, economists etc. It 
can also be used in a statistical sense, such that the most likely events tend to be associated 
with slowly- varying phase variation in the frequency domain, and unlikely events tend to be 
associated with rapidly varying phase with frequency. 

In analyzing mathematical problems in physics, one often finds it desirable to know the 
behavior of a function for large values of some parameter s, that is, the asymptotic behavior 
of the function. Specific examples are furnished by the well-known Gamma function and 
various Bessel functions [42] . All these analytic functions are defined by integrals 



where H is analytic in z and depends on a real parameter s. In the most cases of physical 
interest the integrals like I{s) can be redefined in terms of 
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(1) 



H{z,s) = h{z) exp[sf{z)], 



(2) 
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which leads to the exphcit resolution of I{s), for large values of s, by an asymptotic approx- 
imation procedure named steepest descent method (SDM). The key idea in this method is 
that as s grows very large, the main contribution to the integral will come from values of z 
very close to the points at which df{z)/dz = 0. Irrespective of the nature of the particular 
function, for some paths of integration along which the imaginary part of f{z) is constant, 
these points are called saddle-points because they are the points where the real and imaginary 
parts of f{z) are stationary with respect to the position in the complex plane {x, y) without 
being an absolute maximum or minimum. In fact, if we have f{z) = u{x,y) + iv{x,y), 
we can easily show from Cauchy-Riemann relations that along any path v{x,y) = constant 
the rate of change of u{x,y) can only vanish at a saddle-point [42]. From any point, the 
direction in which u{x, y) decreases most rapidly is one along which v{x, y) is constant, and 
in this sense such paths are paths of steepest descent. Perhaps the most common situation 
where the SDM is applicable is that in which the path of integration v{x,y) = constant 
runs from a saddle-point Zg to infinite, with u{x, y) decreasing monotonically all the way. In 
a simplified analysis, the function f{z) can be expanded around Zg into a Taylor expansion. 
If we replace f{z) with the first two terms of this expansion in the exponential of Eq. (2), 

f{z)^f{zs)-ii\z), (3) 

with 

= Y —{z - Zs), 

and extend the limits of integration from — oo to +00, we find 

/+00 
d^ihizi^i)) exp[-s^i\z)], (4) 
-ex. 

In extensions of this method complex analysis is used to find a contour of steepest descent for 
an equivalent integral expressed as a path integral. To get the first term in the asymptotic 
expansion of (4) the value of fj,{z) h{z)/ f'{z) at z = Zg is required. Provided f"{zs) is not 
zero and h{zs) is not infinite the integral /(s) results 

/ 27r 

H-s) ^ d-jjji^hizs) ey^p[sf{zs)], (5) 

which, however, is not completely free of some particular supplementary conditions for its 
applicabihty [42]. An alternative, and very similar, method to the SDM is the method 
of stationary phase (SPM). Though perhaps less general and less immediately convincing 
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analytically, it often has the advantage of closer contact with the physical problem. The 
integral to be considered is more suitably written as 



where, in practice, the exponential commonly represents a traveling wave. The SPM was first 
introduced explicitly by Lord Kelvin [4]. A rigorous mathematical treatment which would 
justify the statements made above was subsequently done by Watson [43] and the complete 
discussion involving the asymptotic approximation was done by Focke [44]. The deductions 
from the SPM follow much the same pattern as those from the SDM. The assumed paths of 
integration are v{x, y) = constant, which means that the amplitude part of the exponential is 
constant along the path, while the phase part varies most rapidly, a reversal of the situation 
in the SDM. For asymptotic values of s the rapid oscillations of sf{z) over most of the 
range of integration means that the integrand averages to almost zero. Exceptions to this 
cancellation rule occur only when the stationary phase condition is satisfied, i. e. the only 
significant contribution to the integral arises from portions of the path in the vicinity of 
saddle-points or end-points, but the physical interpretation of the mechanism by which this 
comes about is now in terms of phase interference rather than amplitude decay. Thus the 
approximation corresponding to (5) for the SDM is 



for the SPM. But one aspect in which there is some distinction between the methods should 
be noted. With a steepest descent path which starts at a saddle-point and does not go 
to infinity, the contribution of the point at the end of the path to the strict asymptotic 
expansion is zero in comparison with the saddle-point contribution by the virtue of the 
extra exponential factor it contains. On the other hand, with a stationary phase path of the 
same type the contribution of the point at the end of the path is, in general, of the order of 
that of the saddle-point merely divided by y^; it is excluded, therefore, from the asymptotic 
approximation only if the first term of the asymptotic expansion is alone retained. To 
sum up, the methods of steepest descent and stationary phase, strip off their mathematical 
expression, depend on choosing a path of integration in such a way that the integrand, due 
to its exponential factor, contributes irrelevantly to the integral except in the vicinity of 
certain saddle-points (or end-points). 




c 



(6) 
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To illustrate the applicability of the SPM into a physical problem it should be instructive 
to assume an equivalence between the complex variable z and a real one k so that an integral 
which represents a free quantum particle propagation can be identified by the wave packet 
solution of the one-dimensional Schrodinger equation, 



in—tlj{x,t) = 



2m dx^ 

with V{x) — 0, described in terms of the integral, 

dk 



+ V{x) 



tp{x,t) 



G{k, ko) exp [—1 {E{k) t + kx — k Xq)], 



(8) 



(9) 



where we have set h — 1 and the dispersion relation as E{k) — k'^/{2m). The function 
G{k,ko) physically represents a narrow momentum distribution centered around the mo- 
mentum ko, and the Eq. (9) can be identified as the integral of Eq. (6) by simply rewriting 
G{k,ko) as \/2TTh{z) and sf{z) as the phase i{kx — kxo — Et) in an extension to the 
complex plane {k ^ z). 

The integral (9) can therefore be estimated by finding the value for which the phase has 
a vanishing derivative, evaluating (approximately) the integral in the neighborhood of this 
point. The movement of the peak coordinate of the wave packet ip{x, t) can be obtained by 
imposing the stationary phase condition 



-^[Et- k{x-Xo)] 
ak 







Jxig 

m 



(10) 



when ks — ko, the maximum of G{k,ko), which means that the peak of the wave packet 
propagates with a velocity v — In fact, we ratify this result when we explicitly calculate 
the integral of Eq. (9) by introducing a Gaussian momentum distribution given by 

a^ik-koY' 



Gik,ko)=gik-ko)={^'' 



exp 



(11) 



In this case, the result of the integration (9) gives us 



ipix,t) = ip[x - Xo,t], 



(12) 



where 



ip[x,t] = 



2 



exp 



fco ^ 
m 



a' (l 



2it 



arctan (^) -|- Ok^x - E^t) 



which evidently confirms the result of Eq. (10). 



(13) 
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Meanwhile, the method leads to some new interpretive discussion when the momen- 
tum distribution becomes a complex function, i. e. when G{k, ko) can be written as 
|G'(A;, A;„)| exp [i \{k)]. In this case, the stationary phase condition for a free wave packet 
represented by Eq. (10) will be modified by the presence of an additional phase X{k). By 
assuming that X{k) can be expanded in the neighborhood of ko, i. e. 

dX{k) 



X{k) ^ Xiko) + ik- ko) 



dk 

the new stationary phase condition can hereby lead to 



(14) 

k=ko 



d\{k) 



+ -t, (15) 

k=ko 



■~ max ^ u 

which can be interpreted as space-time translation. As it is commonly presented in several 
textbooks of quantum mechanics [45], one can be persuaded to take the stationary phase 
condition as a necessary and sufficient statement for applying the method. Nevertheless, 
although it can be used to predict the most likely outcome of an experiment, even when it is 
related to the applicability of the stationary phase principle, it does not exclude the possibil- 
ity of alternative outcomes. It is well-known in the literature that both SDM and SPM are 
applicable only if certain conditions apply. As we shall demonstrate in the following section, 
the results of the SPM, when applied to the one-dimensional scattering problem, depend 
critically upon the manipulation of the generic amplitude G{k, ko) prior to the application 
of the method. 



A. The above-barrier diffusion problem 

To elucidate some questions upon the ambiguities which appear when we utilize the 

SPM, let us study the scattering of an incoming wave packet by a potential barrier for 
propagating energies larger than the barrier upper limit Vq. The stationary wave solution of 
the Schrodinger equation (8) [78] obtained when we consider the potential barrier described 

by 

if X < region I, 

Vo if < X < L region II, (16) 

ii X > L region 111, 
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can be decomposed into different wave functions for each interval of x, i. e. 



+ 



(f>iiik,x) 



+ 



(pinc{k,x) + (pR{k,x) (t)a{k,x) + 4>D{k,x) 



(j)ui{k,x), 
(t)T{k,x) 



(17) 



where 



4>inc{k,x) = exp[ikx], 

(f)u{k,x) = R{k) exp [—i k x], 

4>a{k,x) = a{k)ex.Y)[iqx\, 

4>f3{k,x) = I3{k)exp[—iqx], 

(j)T{k,x) = T{k)exp[ikx], (18) 

with q — {k^ — w^) 5 and w — \/2m.Vo. By solving the constraint equations [45] which set 
the continuity conditions at x = and x = L, we obtain 



a{k) 
R{k) 



k{k + q) 
Hk) 



Hk) 



exp[iQ{k)-iqL], f3{k) = - 
sin[gL]exp[iG(fe)], T{k) = 



k{k — q) 



Hk) 
2kq' 



exp [i @{k) +iqL], 
exp[ie{k) -ikL], (19) 



where 



and 



J^(fc) = {A k' q^ cos^q L] + (k^ +qy sin' [qL]y 



k' -\- q' 

Q{k) = arctan — — — tan [qL] ^ . 

K q 



(20) 



(21) 



The explicit expression for the correspondent propagating wave packets can be obtained by 
solving integrals hke 

dk 



^f{x,t) 



g{k — ko) (f)f{k,x) exp [—iEt], 



(22) 



with f = a, P, R, T. As a first approximation which is commonly used in quantum me- 
chanics textbooks [45] , we obtain the analytical formulas to these integrals by assuming the 
momentum distribution g{k — ko) is sufficiently sharply peaked around the maximum point 
kg, with k^^ > w. In this case, the integration can be extended from [w, oo] to [— oo, oo] 
without modifying the final result. However, a careful investigation of Q{k) in Eq. (21) 
would clearly indicate that the above integral is not dominated by the exponential for cer- 
tain values of kg/w (see the Fig. (1 where we describe the derivative of Q{k) as a function 
of ko/w). 
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In the sense we are investigating, as we shall demonstrate in the following, it is erroneously 
assumed that the /c-dependent phase terms Q{k) and q = q{k) can be approximated by a 
series expansion around k = up to the first order term, i. e. 

dq 



@{k) ^ Q{ko) + {k- ko)0'{ko) and g go + - ^o) 



dk 



(23) 



k=ko 



where 



e'{k) = ^ 

Q 



{k^ + q^)eqL- {k^ - q^Y sin [q L] cos [q L] 



dq 
dk 



^0 

k=ko ^° 



4fc2 g2 + (^2 _ g2)2 gj^2 [g^] 

At the same time, when k is approximated by ko, it is assumed that the modulating amplitude 
|0/(/co, x)\ can be put out of the integral giving the following results. 



i^Inc{x,t) 



ip[x-Xo,t], 



dk 



g{k — ko) exp [—iEt — i k{x + xq) +i{k — A:o)6'(fco)] 



ipcix,t) w a(fco) 



R{ko) exp [-i ko <d'{ko)] f [-x - + &{ko),t], 
dk 



i^fi{x,t) 



J - fco) exp -i{Et + kxo-qox)+i{k-ko)(^^{x-L) + e'{ko)) 

Q;(fco)exp iqoX - {x - L) - iko@'{ko) tp ^{x - L) - Xq + 0'{ko),t 
/3(fco) j -^g{k-ko) exp[-i{Et + kxo + qox)-i{k-ko)(^^{x-L)-e'{ko)) 



P{ko)exp -iqox + i {x - L) - ikoQ'{ko) 



{x- L)-Xo + e'{ko),t 



Mx,t) « T{ko) 



dk 



g{k — ko) exp [—iEt + i k{x — Xo) +i{k — ko) (©'(fco) — L)] 



T{ko) exp [-i koQ'iko)] (p[x-Xo-L + e'{ko),t]. 



(24) 



By observing the Gaussian shape of ip[x, t] given by Eq. (13), one could easily (but wrongly!) 
identify the position of the peak of the wave packets. The times corresponding to the position 
X of the incident and refiected wave packet peaks would be respectively given by 

x + xo- @'{ky 



tlnc{x) 



X Xn 



Vk 



and tji{x) 



k=ko 



Vk 



(25) 



k=ko 



where Vk = ^ = Only a; < is physical in this result since these waves, by definition, lie 
in region I. Since the phase of the incoming wave contains only the plane wave factors, i.e. it 
is devoid of Q{k), the incoming peak reaches the barrier at x = at time tinc{0) — —{xo/vko) 
(neglecting interference effects). The presence of the phase term Q'{ko) — {dQ{k)/dk)\k=ko 
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for tii{x) would imply a time delay, or a time advance, which depends on the sign of Q'{ko), 
for the reflected wave with respect to the incident one, i.e. 

e'{k) 



tn{0) = tinciO) 



(26) 



k=kQ 



It is analogous to what happens for the step potential tunneling penetration when E <Vq 
[45]. When < x < L in region II, the times for the "intermediary" a and (3 wave packet 
peaks would be given by 



{x - L) [x, - &{k)) 



Vk 



k=ki) 



{x - L) (xo - G'(fc)) 



Vk 



(27) 



and ts{x) = 

And finally, when a; > L in region III, the peak of the transmitted wave packet would be 
written as 



trix) = 



X- xo- L + e'{k) 



(28) 



k=ko 



The above results (25-28) could be directly obtained by applying the SPM to the wave 
functions ipf^Xjt) expressed in Eq. (22), where the stationary phases for each decomposed 
wave component / would be given by 

'&incix,t,k) = -Et + k{x-Xo), 

■dn{x,t,k) = -Et - k{x + xo) + Qik), 

■d^{x,t,k) = -Et-kxo + qix- L) + @{k), 

^is{x,t,k) = -Et - kxo - q{x - L) + @{k), 

'&T{x,t,k) = -Et + k{x-Xo-L) + @{k). (29) 

Meanwhile, these time-dependencies must be carefully interpreted. In the Fig. (2), we il- 
lustrate the wave packet diffusion analytically described by the results of Eq. (24). The 
"pictures" display the wave function in the proximity of the barrier for suitably chosen 
times. By taking separately the right {a) and left (/?) moving components in region II, 
independent of the value of Q{k), we can observe that 



L 



(30) 



k=ko 



Ai„ = t„{L) - t„(0) = At0 = t^iO) - tp{L) = 

which correspond to transit time values "classically" expected. 

However, by considering the information carried by the wave packet peaks, we would 
have a complete time discontinuity at x — represented by Eq. (26) and by the fact that 



Xo ^ -L _ e^(fc) 

Vk Vq Vk 



/ tlnc{0), 



(31) 



k=ko 
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i. e. the a wave could appear in region II at a time t„(0) before (or after) the incident 
wave arrives at x = 0. The results expressed by Eqs. (26) and (31) which are illustrated in 
Fig. 2 immediately imply the non-conservation of probabilities since they lead to disconti- 
nuity points (x = and/or x = L) for the wave functions, which are Schrodinger equation 
solutions, , and for its spatial derivative. The normalized squared modulus of this wave 
functions represents the spatial probability distribution of finding the corresponding propa- 
gating particle. At x = (and/or x — L) the incoming flux of probabilities is different from 
the outgoing flux, which, from basic definitions of quantum mechanics continuity equation 
[45] , leads to the non-conservation of probabilities. It is, by principle, unacceptable. 

To be more accurate with such analysis and clear up some dubious points, it should be 
more convenient to verify two simple particular cases. In order to simplify the calculations, 
we set Xo — and we choose QoL — mr {n — 1, 2, 3, ...) so that we can write, 

e(A;o) = 0, Riko) = 0, T{ko) = l (32) 

and 

e\k) = !^tl±A > aili (33) 

Qo 2koqo Qo 



so that 



^2 



U0) = ^%^>0. (34) 



Otherwise, if we choose QqL = {n + |)7r, we shall obtain, 
and 

= rnL ^ ruL 

^ ' qo kl + ql qo ^ ' 

which gives 

mL(|-*£^^ (37) 

The latter negative vahic ilhistrativcly corroborates with the absurd possibility of the ap- 
pearance of peak associated with the a wave occurring before the arrival of the peak of the 
incident wave packet at x = 0, a clear representation of discontinuity at x = 0. As we have 
stated before, the phase derivative Q'{k) illustrated in Fig.(l) does not present the adequate 
behavior for applying the SPM accurately. In spite of this, it is currently ignored in the 
literature. 
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B. Analytic solution for multiple peak decomposition 



Differently from the analytical analysis presented above, the numerical simulations of a 
wave packet diffusion by the potential barrier shows the appearance of multiple peaks due 
to the two reflective interfaces at x = and x = L (see the figures). In fact, numerical 
calculations automatically preserve probabilities, at least to within the numerical errors. 
This observation suggest a new analysis and subsequent interpretation/quantification of the 
ambiguities presented in the previous section. In order to correctly apply the SPM and 
accurately obtain the position of the peak of the propagating wave packet with an accurate 
time dependence, we are constrict to solve the continuity constraints of the Schrodinger 
equation at each potential discontinuity point a; = and x = L hj considering multiple 
successive reflections and transmissions. The phenomenon was already described in [16] and 
consists in an incoming wave of unitary amplitude with momentum distribution centered 
at ko which reaches the interface x — where, at an instant t — to, is decomposed into 
a reflected wave component of amplitude (see the diagram) and a transmitted wave 
component of amplitude a^. The transmitted wave continues propagating until it reaches 
X — L where, at t — t^, it is now decomposed into a reflected wave component of amplitude 
/3i and a transmitted wave component of amplitude T^. That new reflected wave component 
comes back to x — where, at t — t2, it is again decomposed into another reflected wave 
component of amplitude 0:2 and transmitted wave component of amplitude R2 which will be 
added to R^. This iterative process continues for infinity times as we can see in the following 
diagram. 



exp[i kx] + Ri cxp[—i k x] 
i?2 exp[— i kx] 

R„ exp[—ikx] 
x = 



a2exp[iqx] + I3iexp[—iqx] 

: + : 
a„ exp[i qx]+ /3„_i exp[-z q x] 



Q;icxp[iga;] + l3iexp[—i q x] 
a2exp[iqx] + P2exp[—iqx] 

: + : 
Q!„ exp[i qx]+ Pri exp[-i q x] 
X = L 



Ti explikx] 
T2 exp[ikx] 

T^explikx] 
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The continuity constraints over ip{x, t) for each potential step at a; = and x — L determine 
the coefficients 



k — q 

ikq 



2k ^ 2k{q-k) 



(/j + ®^P[^ ~ ^) ^] ' R2 = -j^Oil(3l 



and estabhsh the recurrence relations 



R, 



+2 



k — q 



exp[2iq L] 



n = l,2,... 



Rn+l /?n + Q y 

which allow us to write the sum of the coefficients a„,(3„, and as 

' k — q^ ^ 



(38) 



(39) 



R 



a 



R„ — Rl + R; 

n=i 

oo 

y ^ an = ai 

n=i 

oo 

n=i 

oo 



1 



exp[2i(/L] 



k — q 
k + q 



k + q 

k — q^ ^ 
k + q 



exp[2igL] 
exp[2igL] 
exp[2 iqL 



(40) 



The above summations reproduce exactly the expressions in Eq. (19). In this form the 
interpretation is easy represents the first reflected wave (it has no time delay since it 
is real). R2 represents the second reflected wave and, as a consequence of the continuity 
condition at a; = 0, it is the sum, in region II, of the first left-going wave (/3i) and the second 
right-going amplitude (q;2), i.e. 

q 

R2 = a2 + Pi = — ai l3i . 

k 

This structure is that given by considering two "step functions" back-to-back. Thus at each 
interface the "reflected" and "transmitted" waves are instantaneous i. e. without any delay 
time. 

Now we can calculate the time at which each transmitted and/or reflected wave appear by 
applying the SPM for each component of the total transmitted T or reflected R coefficients. 
It will give us the recurrence relation 

m{x — L) 



tn — ^n— 1 "I" 



(41) 



which coincides with the "classically" predicted value for the velocity of the particle above 
the barrier. Indeed the SPM applied separately to each term in the above series expansion 
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for R yields delay multiples of 2 {dq/dE)\^^^^ L — 2{rn/qo) L. This agrees perfectly with 
the fact that since the peak momentum in region II is q^, the a and (5 waves have group 
velocities of qo/m and hence transit times (one way) of (m/qo) L. The first transmitted peak 
appears (according to this version of the SPM) after a time {m/qo) L, in perfect agreement 
with the above interpretation. 

Is this compatible with probability conservation? It is because of the following identity 

oo 

^i\R„\' + \T„\') = l . (42) 

n=i 

This result is by no means obvious since it coexists with the well-known result, from the 
plane wave analysis, 

oo oo 

\R\' + \T\' = \^Rn\' + \^T„\' = l. (43) 

n=i n=i 

To conclude, we can state that the conditions for applying the SPM in scattering problems 
of such a type depend upon the correct manipulation of the amplitude. A posteriori this 
seems obvious but the point is that unless we know the number of separate peaks the SPM 
is inaccurate. There is, however, a converse to this question. If the modulating function is 
such that two or more wave packets overlap, then we cannot treat them separately without 
considering all the interference effects. The above-barrier analysis is simply a particular 
example of this ambiguity. 

C. The analytical formula 

Using the linear approximation (23) for the momentum q and the Gaussian momentum 
distribution (11), and considering the expressions in (40) resulting from the multiple peak 
decomposition in (22), newly with the same pertinent approximations used for obtaining 
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(24), we can analytically construct the following simplified expressions 



ipi{x,t) 

tljR{x,t) 
ll)c,{x,t) 

i>T{x,t) 



ip[x - Xo, t], 

oo 
n=0 

oo 



n=0 



2fco 

ko+qo 



90 



— 2fco(TO-fco) 



go 



oo 

( |o_2o exp [-2 l' 

\*:o+9() I <?() 

n=() 

oo 



) ^[-a;-a;o + 2(n+l)^L, t], 



, 2n 



ip[{x + 2nL)^-Xo,t], 



OO 

n=0 



10 



)2n 
ip[{2nL 



+ 2L-x)!^-Xo, t], 



OO 

-i 2"! l] y f exp 

n=0 



QO 



)2n 
ip[x- 



x,-L+{2n+l)^L 



,,t]. 



(44) 



which allow us to visualize the multiple peak decomposition illustrated in Fig. 3. 

Comparing with the exact numerical results, the analytic approximations expressed in 
(44) will be valid only under certain restrictive conditions. To explain such a statement, 
when we analytically solve the integrals like (22), the exponential exp [2iqL] which appears 
in the recurrence relations (39) are approximated by 



k 

exp [2i Qo L + 2i {k — ko)—]. 



(45) 



In this case, the above analytical expressions can be obtained only when the exponential 
function does not oscillate in the interval of relevance of the envelop Gaussian function 
g{k — ko), i. e. when AqL < tt, which can be expressed in terms of the wave packet width 
as 

ko L 



AqL 



dq 
dk 



AkL K — - < TT. (46) 
k=ko 1o « 

where we have used Ak ~ 1/a. The above constraint implicitly carries a very peculiar 
character: the multiple peak decomposition can also be evidenced when the wave packet 
width a is larger than the potential barrier width L, i. e. when we assume the relation ko/qo 
is in the interval 1 < ko/qo < vr we can have propagating wave packets with a > L satisfying 
the requirements for the multiple peak resolution (44). To illustrate this possibility, in Fig. 4 
we plot the density of probabilities representing the collision of a wave packet of average 
width a with a potential barrier Vq of width L — 0.8 a. We compare the analytical (line) 
with the numerical (star symbols) results from which we can evidently observe the accurate 
equivalence between them since we are respecting the condition (46). Both the figures (3) 
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and (4) display the square of the wave function in the proximity of the barrier for suitably 
chosen times. One clearly sees the appearance of multiple peaks due to the two reflective 
interfaces at x — and x — L. 



III. LIMITATIONS ON THE TUNNELING ANALYSIS 



Generically speaking, the SPM can be successfully utilized for describing the motion of 
the center of a wave packet constructed in terms of a momentum distribution g{k — ko) which 
has a pronounced peak around kg. In case of a potential barrier given by, 



V{x) 



X £ [-L/2, L/2] 



(47) 



[-L/2, L/2] 

it is well known that the transmitted wave packet solution {x > L/2) calculated by means 
of the Schrodinger formalism is given by [45] 



tp^ix, 



dk 



g{k-ko)\T{k,L)\ exp 



where, in case of tunneling, the transmitted amplitude is written as 



ik{x- L/2) -i-z—t + i @{k, L) 



|T(fc,L)|= 1 + 



w 



4A;2 

and the phase shift is described in terms of 



sinh^ [p{k) L] 



Q{k, L) = arctan 



2e 



w 



kp{k) 



tanh[p(A;)L] L 



(48) 



(49) 



(50) 



for which we have made explicit the dependence on the barrier length L and we have adopted 
p{k) = (w^ — /c^)^ with w = (2mVo)2 and h = 1. Without thinking over an eventual 
distortion that \T{k,L)\ causes to the supposedly symmetric function g{k — ko), when the 
stationary phase condition is applied to the phase of Eq. (49) we obtain 



x-L/2 



k„ 



dQ{k,L) 



0. 



(51) 



m dk 

The above result is frequently adopted for calculating the transit time tx of a transmitted 
wave packet when its peak emerges at x — L/2, 



m dG(A;, a(i)) 



dk 



2m L \ w'^ sinh {a) cosh (a) — (2 A;- 



,2 

max 



4 {w^ - klJ + sinh^ (a) 



(52) 
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where we have defined the parameter a — {w"^ — k^^^)^ L. The concept of opaque hmit is 
introduced when we assume that k^^^ is independent of L and then we make a tends to oo 
[11]. In this case, the transit time can be rewritten as 

OL 2m 

^ — ^(T^- ^^^^ 

"-max PV'^max / 

In the literature, the value of k^^ is frequently approximated by ko, the maximum of g{k—ko), 
which, in fact, does not depend on L and could lead us to the transmission time superluminal 
interpretation [11, 46, 47] of (53). To clear up this point, we notice that when we take the 
opaque limit with L going to oo and w fixed as well as with w going to oo and L fixed, with 
ko < -u; in both cases, the expression (53) suggests times corresponding to a transmission 
process performed with velocities higher than c [11]. Such a superluminal interpretation was 
extended to the study of quantum tunneling through two successive barriers separated by a 
free region where it was theoretically verified that the total traversal time does not depend 
not only on the barrier widths, but also on the free region extension [47]. Besides, in a 
subsequent analysis, the same technique was applied to a problem with multiple successive 
barrier where the tunneling process was designated as a highly non-local phenomenon [46]. 

It would be perfectly acceptable to consider k^^^ = k„ for the application of the stationary 
phase condition when the momentum distribution g{k — ko) centered at ko is not modified by 
any bound condition. This is the case of the incident wave packet before colliding with the 
potential barrier. Otherwise, the procedure for obtaining all the above results for the the 
transmitted wave packet does not take into account the bounds and enhancements imposed 
by the analytical form of the transmission coefficient. 

Hartman himself had already noticed that the 

" Discussion of the time dependent behavior of a wave packet is complicated by its dif- 
fuse or spreading nature; however, the position of the peak of a symmetrical packet can be 

described with some precision" 

The origin of the Hartman effect, and the above- illustrated saturation of group delay with 
barrier length, has been a mystery for decades. The phenomenon exists in the tunneling of 
all kinds of waves, be they matter waves, electromagnetic waves, or sound waves. Experi- 
mentally it has been observed using electromagnetic waves and sound waves and there is no 
question as to its existence. If one interprets the group delay in tunneling as a transit time 
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then the Hartman effect naturally leads to superluminal and unbounded group velocities. 
However, to perform the correct analysis, we should calculate the right value of k^^^ to be 
substituted in Eq. (52) before taking the opaque limit. We have to consider the relevant 
amplitude for the transmitted wave as the product of a symmetric momentum distribution 
g{k — ko) which describes the incoming wave packet by the modulus of the transmission am- 
plitude T{k, L) which is a crescent function of k. The maximum of this product representing 
the transmission modulating function would be given by the solution of the equation 



g{k-K) \T{k,L)\ 



g'{k-ko) ^ |r(A;,L)|' 



= (54) 



[g{k-ko) \T{k,L)\ 

Obviously, the peak of the modified momentum distribution is shifted to the right of so 
that /Cmax have to be found in the interval ]ko,w[. Moreover, we confirm that k^^ presents 
an implicit dependence on L as we can demonstrate by the numerical results presented in 
Table 1 where we have found the maximum of g{k — ko) \T{k, L)\ by assuming g{k — ko) is 
a Gaussian function almost completely comprised in the interval [0, w] given by 



g{k -ko) = (^] exp 



,2 



4 



a\k-ko)^^ 



(55) 



2tt^ 

By increasing the value of L with respect to a, the value of /c^ax obtained from the numerical 
calculations to be substituted in Eq. (52) also increases until L reaches certain values for 
which the modified momentum distribution becomes unavoidably distorted. In this case, 
the relevant values of k are concentrated around the upper bound value w. We shall show 
in the following that the value of L which starts to distort the momentum distribution can 
be analytically obtained in terms of a. 

Now, if we take the opaque limit of a by fixing L and increasing w, the above results 
immediately ruin the superluminal interpretation upon the result of Eq. (52) since t^^ tends 
to oo when k^^^ is substituted by w. 

Otherwise, when w is fixed and L tends to oo, the parameter a calculated at k — w 

becomes indeterminate. The transit time tr still tends to oo but now it exhibits a peculiar 

dependence on L which can be easily observed by defining the auxiliary function 

sinh (a) cosh (a) - a 

G{a) = ■ ,2/ X (56) 

Sinn (aj 



which allow us to write 



t^ = ^JI^G{a). (57) 

wa 
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When a ^ 1, the transmission time always assume infinite values with an asymptotic 
dependence on {w^ — k^)~^ , 

oo. (58) 



w {w'^ — A;2) 2 

Only when a tends to we have an explicit linear dependence on L given by 



Q '2 m L ( G{a)\ AmL 



t^ = ^lim^ =-— (59) 

In addition to the above results, the transmitted wave must be carefully studied in terms 
of the rapport between the barrier extension L and the wave packet width a. For very thin 
barriers, i. e. when L is much smaller than a, the modified transmitted wave packet presents 
substantially the same form of the incident one. For thicker barriers, but yet with L < a, 
the peak of the Gaussian wave packet modulated by the transmission coefficient is shifted 
to higher energy values, i. e. k^^^ > ko increases with L. For very thick barriers, i. e. 
when L > a, we are able to observe that the form of the transmitted wave packet is badly 
distorted with the greatest contribution coming from the Fourier components corresponding 
to the energy w just above the top of the barrier in a kind of filter effect. We observe that 
the quoted distortion starts to appear when the modulated momentum distribution presents 
a local maximal point at k — w which occurs when 

d 



[g{k-K) \T{k,L)\\ 



> 0. (60) 

k=w 



Since the derivative of the Gaussian function g{k — ko) is negative at k — w, the Eq. (60) 
gives the relation 



^'^^ - < lim 



T'ik,L) 



g{w — ko) k^w [ r(fc, L) _ 
which effectively represents the inequality 



(l + ^) 



< ^ (61) 



/ , N wL"^ ^ /3 / ko\ 

-in:-ko)<^ L>^-a[l-^y (62) 

Due to the filter effect, the amplitude of the transmitted wave is essentially composed 
by the plane wave components of the front tail of the incoming wave packet which reaches 
the first barrier interface before the peak arrival. Meanwhile, only whether we had cut 
the momentum distribution off at a value of k smaller than w, i. e. /c ~ (1 — 5)w, the 
superluminal interpretation of the transition time (53) could be recovered. In this case, 
independently of the way as a tends to oo, the value assumed by the transit time would 
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be approximated by tj^ ~ 2m/ w 5 which is a finite quantity. Such a finite value would 
confirm the hypothesis of superluminality. However, the cut off at k (1 — 5)w increases 
the amplitude of the tail of the incident wave as we can observe in Fig. 5. It contributes so 
relevantly as the peak of the incident wave to the final composition of the transmitted wave 
and creates an ambiguity in the definition of the arrival time. 

To summarize, at this point we are particularly convinced that the use of a step- 
discontinuity to analyze signal transmissions in tunneling processes deserves a more careful 
analysis than the immediate application of the stationary phase method since we cannot find 
an analj^ic-continuation between the above barrier case solutions and the below barrier case 
solutions. A suggestive possibility can ask for the use of the multiple peak decomposition 
technique developed for the above-barrier diffusion problem [16] presented in section III. 
Thus, in the same framework, we suggest in the following section a suitable way for com- 
prehending the conservation of probabilities for a very particular tunneling configuration 
where the asymmetry presented here is eliminated, and the phase times can be accurately 
calculated. 

IV. THE SYMMETRICAL TUNNELING CONFIGURATION 

Due to the factual influence of the amplitude of the transmitted wave in determining 
the right position of a tunneling particle, we can introduce an alternative analysis where 
the possibility of multiple peak decomposition [16] may be taken into account. By means 
of such an experimentally verifiable exercise, we shall be able to understand how the filter 
effect can analj^ically affect the calculations of transit times in the tunneling process. 

In order to recover the scattered momentum distribution symmetry conditions for ac- 
curately applying the SPM, we assume a symmetrical colliding configuration of two wave 
packets traveling in opposite directions. By considering the same barrier represented in 
(47) , we solve the Schrodinger equation for a plane wave component of momentum k for two 
identical wave packets symmetrically separated from the origin a; = 0, at time of collision 
t = —{mL)/{2ko) chosen for mathematical convenience, where we assume they perform a to- 
tally symmetric simultaneous collision with the potential barrier. The wave packet reaching 
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the left(right)-side of the barrier is represented by 

/ + 00 
dk g{k - K)<t)''^''^ {k, x) exp [-i E t] (63) 
-oo 

where we assume the integral can be naturally extended from the interval [0, w] to the 
interval [—00, +00] as a first approximation. In the framework of the multiple peak de- 
composition [16], we have suggested a suitable way for comprehending the conservation of 
probabilities where the asymmetric aspects previously discussed [65] could be totally elimi- 
nated. By considering the same rectangular barrier V{x), we solve the Schrodinger equation 
for a plane wave component of momentum k for two identical wave packets symmetrically 
separated from the origin a; = 0. By assuming that (^^^"^k, x) are Schrodinger equation 
stationary wave solutions, when the wave packet peaks simultaneously reach the barrier (at 
the mathematically convenient time t — — (mL)/(2/co)) we can write 

^f(«)(fc,x) = exp[±i/ca;] -f-i?^(«)(/c,L)exp[:^i/cx] x < -L/2{x > L/2), 

(/)^(^)(A:,x) = < (j)';^''\k,x) = 7^(«)(fe)exp [^px] + /?^(«)(fe) exp [±px] - L/2 < x < L/2, 

(^3^''^(A;,x) = r^(«)(/c,L)exp[±iA;x] x > L/2 (x < -L/2). 

where the upper(lower) sign is related to the index L{R) corresponding to the incidence on 
the left(right)-hand side of the barrier. By assuming the conditions for the continuity of 0^'^ 
and their derivatives at a; = —L/2 and x = L/2, after some mathematical manipulations, 
we can easily obtain 



^' ^ ^ \l-exp[2/9(A;)L]exp[i2 6l(A;)] J 



(64) 



and 



and R^^'^^k, L) and T'^^^\k, L) are intersecting each other. 



where 



A. Symmetrized and antisymmetrized quantum tunneling configuration 



Since the above introduced colliding configuration is spatially symmetric, the symmetry 
operation corresponding to the 1 <-> 2 particle exchange can be parameterized by the position 
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coordinate transformation x — > —x. At the same time, it is easy to observe that 

</)^(«)(fc,x) = ^[it^{k,x) = 4>itz{K-x) = </)«(^)(A;,-x) (67) 

where the L ^ R interchange is exphcit. Consequently, in case of analyzing the collision of 
two identical bosons, we have to consider a symmetrized superposition of the L and R wave 
functions, 

(^+(A;, x) = (p\k, x) + (p^ik, x) = (l)''{k, -x) + (t)\k, -x) = 4)+{k, -x). (68) 

Analogously, in case of analyzing the collision of two identical fermions (just taking into 

account the spatial part of the wave function), we have to consider an antisymmetrized 
superposition of the L and R wave functions, 

<l)-{k,x) = (f)''{k,x) - 4>''{k,x) = 4>''{k,-x) - <^^(A;,-x) = -4>_{k,-x). (69) 

Thus the amplitude of the re-composed transmitted plus reflected waves would be given by 
R'"'"'{k, L)-\-T'^'^{k, L) for the symmetrized wave function and by R'"'^{k, L)—T^''"{k, L) for 
the antisymmetrized wave function 0_. Reporting to the multiple peak decomposition [16] 
applied to such a pictorial symmetrical tunneling configuration [58], we can superpose the 
amplitudes of the intersecting probability distributions before taking their squared modulus 
in order to obtain 

= exp{-i[/cL-(^±(/c,L)]} (70) 

with 

2fcp(fc) smh[/>(fc)L] ) 
{K' ~ p'{k))cosh[p{k)L\±w'']' ^ ' 

where the plus sign is related to the results obtained for the a symmetrized superposition 
and the minus sign is related to the symmetrized superposition. Independently of the odd 
or even symmetrization of the wave function, we observe from Eq. (70) that, differently 
from the previous standard tunneling analysis, by adding the intersecting amplitudes at 
each side of the barrier, we keep the original momentum distribution undistorted since 
|i?^'^(/c,L) ±T^'^(/c,L)| is equal to one. The previously pointed out incongruities which 
cause the distortion of the momentum distribution g{k — k) are completely eliminated and 
we recover the fundamental condition for the applicabihty of the SPM for accurately de- 
termining the position of the peak of the reconstructed wave packet composed by reflected 
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ip±{k, L) = — arctan 



and transmitted superposing components. The phase time interpretation can be, in this 
case, correctly quantified in terms of the analysis of the novel phase (p±{k, L) since the novel 
scattering amplitudes g{k — ko)\R^'"- ± T'^'^\ ^ g{k — k^) maintains its previous symmetrical 
character. The transmitted and reflected interfering amplitudes results in a unimodular 
function which just modifies the envelop function g{k — k) by an additional phase, and the 
scattering phase time results in 



^(a) m d(p{k, a) 2m L sinh (a) ziz ak 
'^''^'^ k dk ka 2k^ — w'^ ±w'^ cosh (a) 



(72) 



where k — > kg, with a previously defined. The old phase Q{k, L) (Eq. 50) appears when 
we treat separately the momentum amplitudes T{k,L) and R{k,L), which destroys the 
symmetry of the initial momentum distribution g[k — k^) by the presence of the multiplica- 
tive term T{k,L) or R{k,L), and the novel phase ifi{k,L) appears only when we sum the 
tunnneling/scattering amplitudes so that the symmetrical character of the initial momentum 
distribution is recovered (due to the result of Eq. (70). 

We should assume that in the region inside the potential barrier, the refiecting and trans- 
mitting amplitudes should be summed before we compute the phase changes. Obviously, 
it would result in the same phase time expression as represented by (72). In this case, the 
assumption of there (not) existing interference between the momentum amplitudes of the 
refiected and transmitted waves at the discontinuity points x — —L/2 and x — L/2\s purely 
arbitrary. Consequently, it is important to reinforce the argument that such a possibility 
of interference leading to different phase time results is strictly related to the idea of using 
(or not) the multiple peak (de) composition in the region where the potential barrier is lo- 
calized. To illustrate the difference between the standard tunneling phase time and the 
symmetrica/ scattering phase time we introduce the parameter n — k'^ jvf' and we define 
the classical traversal time = {mL)/k. In this case, we can obtain the normalized phase 
times 

(a) 2rfe J cosh (a) sinh (a) — a n (2n — 1) ]^ 

""^l [4n (1 - n) + sinh^ (a)] J ^ ' 

and 

^(a) _ 2rfc [ na± sinh (g) \ ^^^^ 



^''^i a \ 2n - 1 ± cosh (a) 
In order to illustrate the difference between the standard tunneling phase time tT and the 
alternative scattering phase time t'^ we introduce the new parameter n — kl^^/w"^ and the 
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classical traversal time r = (mL) / k^^^ in order to define the rates 

RT{a) = ^ = ^ r cosh(.) sinh(a) --n {2n - 1) j 
^ ' T a\ [4n (1 - n) + sinh^ (a)] J ^ 

and 

2 f n a + sinh (a) 1 
= _L = _ ; k (76) 

^ r a \ 2n - 1 + cosh (a) J ^ 
which are plotted in the Fig. (6) for some discrete values of n varying from 0.1 to 0.9, from 
which we can obtain the most common limits given by 

lim {B^{a)] = lim {i?r(a)} = (77) 

a-^oo a— >oo 

and 

lim {i?T(a)} = 1 + ^, lim {Rl{a)] = 1 + - (78) 

Both of them present the same asymptotic behavior which, at first glance, recover the 
theoretical possibility of a super luminal transmission in the sense that, by now, the SPM can 
be correctly applied since the analytical limitations are accurately observed. At this point, it 
is convenient to notice that the superluminal phenomena, observed in the experiments with 
tunneling photons and evanescent electromagnetic waves [5, 22, 23, 48], has generated a lot 
of discussions on rclativistic causality. In fact, superluminal group velocities in connection 
with quantum (and classical) tunnelings were predicted even on the basis of tunneling time 
definitions more general than the simple Wigner's phase time [19] (Olkhovsky et al., for 
instance, discuss a simple way of understanding the problem [12]). In a causal manner, it 
might consist in explaining the superluminal phenomena during tunneling as simply due to 
a reshaping of the pulse, with attenuation, as already attempted (at the classical limit) [49], 
i. e. the later parts of an incoming pulse are preferentially attenuated, in such a way that 
the outcoming peak appears shifted towards earlier times even if it is nothing but a portion 
of the incident pulse forward tail [5, 50]. In particular, we do not intend to extend on the 
delicate question whether superluminal group-velocities can sometimes imply superluminal 
signalling, a controversial subject which has been extensively explored in the literature about 
the tunneling effect ([12] and references therein). 

Turning back to the scattering time analysis, we can observe an analogy between OTir 
results and the results interpreted from the Hartman Effect (HE) analysis [10]. The HE 
is related to the fact that for opaque potential barriers the mean tunneling time does not 
depend on the barrier width, so that for large barriers the effective tunneling-velocity can 
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become arbitrarily large, where it was found that the tunneling phase time was independent 
of the barrier width. It seems that the penetration time, needed to cross a portion of a 
barrier, in the case of a very long barrier starts to increase again after the plateau corre- 
sponding to infinite speed proportionally to the distance [79]. Our phase time dependence 
on the barrier width is similar to that which leads to Hartman interpretation as we can infer 
from Eqs. (77-78). Only when a tends to we have an explicit linear time-dependence on 
L given by 



which agree with calculations based on the simple phase time analysis where tx — 



compute the phase times illustrated in the Fig. (6) are not effectively constructed with the 
same momentum distributions. The phase Q{k,L) appears when we treat separately the 
momentum amplitudes g{k — k^) \T{k, L)\ and g{k — ko)\R{k, L)\ and the other one ip{k, L) 
appears when we sum the amplitudes of the same side of the barrier. In this sense, as a 
suggestive possibility for partially overcoming the incongruities (here pointed out and quan- 
tified) which appear when we adopt the SPM framework for obtaining tunneling phase times, 
we have claimed for the use of the multiple peak decomposition [16] technique presented in 
the study of the the above-barrier diffusion problem [16]. We have essentially suggested a 
suitable way for comprehending the conservation of probabilities for a very particular tun- 
neling configuration where the asymmetry presented in the previous case was eliminated, 
and the phase times could be accurately calculated. An example for which, we believe, we 
have provided a simple but convincing resolution. 

V. EXACT CORRESPONDENCE BETWEEN PHASE TIMES AND DWELL 
TIMES 

At this point, one could say metaphorically that two bosonic particles represented by the 
symmetrized incident wave function spend a time equal to iT,v3+ inside the barrier before 
retracing its steps or tunneling and that two fermionic particles represented by antisym- 
metrized incited wave function spend a time equal to tr.v?- The answer is in the definition 
of the dwell time for the same colliding configuration which we have proposed. The dwell 




(79) 



2mL 



(l -|- However, it is important to emphasize that the wave packets for which we 
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tD,± = ^ / dx\ct>±,^{k,x)\' (80) 



time is a measure of the time spent by a particle in the barrier region regardless of whether 
it is ultimately transmitted or reflected [26] , 

-L/2 

where ji„ is the flux of incident particles and (l)2{k,x) is the stationary state wave function 
depending on the colliding configuration that we are considering (symmetrical or standard). 
To derive the relation between the dwell time and the phase time, we reproduce the vari- 
ational theorem which yields the sensitivity of the wave function to variations in energy. 
After some elementary manipulations of the Schrodinger equation [51], we can write 

. 1 d /50 50t \ 

Upon integration over the length of the barrier we find 



(82) 

-l/2 \UI1/ UX UHiUX J 

In the barrier limits [x = ±L/2), for the symmetrical configuration that we have proposed, 
we can use the superposition of the scattered waves to substitute in the right-hand side of 
the above equation, 

^ ^[^''\k,x)±ct>f'\k,x) 

L/2(+t/2) ^ 

= ex.Y>[r^ikx\+ eyc£)\^ikx + i{ip±{k,L) — kL)\ (83) 



^±(fc,x)|_ 



By evaluating the right-hand side of the Eq. (83), we obtain 

dk dip+ m /■+^/^ Imlexp (iip+)] dk 

8E^ = T 4^al- 

The first term of the above equation is the phase time or the aforementioned group delay 
i^^. The second term leads to the exphcit computation of the dwell time. By respecting the 
continuity conditions of the Schrodinger equation solutions, in the barrier region we obtain 

a stationary wave symmetrical or antisymmetrical in x, 



= I 7) ^^^^^ ^ ^^'^ ^^^^^ 



(85) 



which, from Eq. (80), leads to 



(q) ^ 2Tfcn Q;zfcsinh(a) 

^■^± a 2n-l±cosh(a) ^ ' 
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The self-interference term which comes from the momentary overlap between the incident 

and the reflected waves in front of the barrier is given by [80] 

.(a) ^ Jm[exp {iip±)] dk ^ m sin ^ 2rfc (l-n)sinh(a) 

~ k dE~ k^ ~ a 2n-l±cosh(a)' ^ ^ 

It is interesting to observe that the result for the self-interference delay t^i^^^ in the above 
equation also depends on the parity of the wave function. The dwell time is obtained from 
a simple subtraction of the quote self-interference delay ti^^,^ from the phase time that 
describes the exact position of the peak of the scattered wave packets, tT,ip± — ^D,(p± ~l~ ^i,>p± 
as we can notice in the Fig. 7. 

Adopting the classical traversal time = {mL) / k for normalizing the results illustrated 
in the above figures allows us to notice an important aspect concerned with the two identical 
fermion collision. The two identical fermion collision leads to the possibility of an unusual 
accelerated tunneling transmission. 

In fact, when each separated transmission coefficient T^'^ prevails over the each reflection 
coefficient i?^'^, i. e. |T|^ > |it!|^, we have |T|^ > 1/2. For satisfying such a requirement 
the Eq. (49), after some manipulations, leads to {wL)/{2y/n) < wLsmh{a))/{2y/n) < 1. 
For two identical bosonic particles, the possibility of accelerated tunneling transitions with 
respect to the traversal classical course is quantified by t'^l,^ < It occurs only when 
a/2 > {a/2) tanh {a/2) > 1. Since a — wL\/l — n, the intersection of the "weak version" of 
both of the above constraints, {wL)/{2^Jn) < 1 and a/2 > 1, leads to n > 2, which definitely 
does not correspond to an effective tunneling configuration. In the region where the one-way 
direction transmission coefficient prevails over the reflection coefficient, bosons should tunnel 
with a retarded velocity with respect to the classical velocity since, in this case, > r^. 
The supposed accelerated transit of the tunneling wave packet, and therefore, superluminal 
velocities and the Hartman effect, would never occur for |T|^ > \RY. It does not correspond 
to the theoretical result for two fermionic particles since we can easily notice that ^t^^. is 
always smaller than r^,. In fact, ^t"^. < t^. does not have a correspondence with the it! and 
T coefficients. Consequently, the possibility of accelerated tunneling and the verifiability of 
the Hartman effect for the two identical fermion tunneling configuration is, in fact, concrete. 

Claiming for the relevance of the use of the multiple peak decomposition technique, we 
have obtained the transit times for a symmetrized (two identical bosons) and an antisym- 
metrized (two identical fermions) quantum tunneling configuration. We have demonstrated 
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that the phase time and the dwell time give connected results in spite of the exact position 
of the scattered particles being explicitly given by the phase time (group delay). For the 
antisymmetrized wave function configuration, an unusual effect of accelerated tunneling ef- 
fect have been clearly identified in simultaneous two identical fermion tunneling. In spite 
of quoting the superluminal interpretation, our discussion concerned with the definition of 
the strict mathematical conditions stringing the applicability of the stationary phase prin- 
ciple in deriving transit times. Even with the introduced modifications, our results partially 
corroborate with the analysis [13, 21] that gives an answer to the paradox of the Hartman 
interpretation [10] since we provide a way of comprehending the conservation of probabilities 
[16, 65] for a very particular tunneling configuration where the asymmetry and the distor- 
tion aspects presented in the standard case were all eliminated. Otherwise, one should keep 
in mind that accelerated tunneling transmission and, generically, the Hartman Effect, even 
in its more sophisticate consequences appears to have been experimentally verified [74], in 
particular, for opaque barriers and nonresonant tunneling [75], and, under severe analytical 
restrictions, reproduced also by numerical simulations and constrained theoretical analysis 
[64, 76]. In this sense, we suggest that the results here obtained deserve further attention by 
experimenter since the existing experiments, under the theoretical focus, report only about 
inaccurate results derived from non-symmetrical quantum tunneling constructions. 



VI. DELAY TIMES FOR RELATIVISTIC TUNNELING 

It is very difficult and probably even confusing to treat all interactions of plane waves 
or wave packets with a barrier potential using a relativistic wave equation [53-56]. This is 
because the physical content depends upon the relation between the barrier height and 
the mass m of the incoming (particle) wave, beside of its total energy E. In the first attempt 
to evaluate this problem, Klein [59] considered the refiection and transmission of electrons 
of incidence energy E on the potential step V{x) — Q{x)Vo in the one-dimensional time- 
independent Dirac equation which can be represented in terms of the usual Pauli matrices 
[60] by [81] 

[a^a'di -{E- @{x)Vq) - a^m] (f){k, x) = 0, (from this point c = h=l), (88) 
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which corresponds to the reduced representation of the usual Pauh-Dirac gamma matrix 
representation. The physical essence of such a theoretical configuration lies in the prediction 
that fermions can pass through large repulsive potentials without exponential damping. It 
corresponds to the so called Klein tunneling phenomenon [54] which follows accompanied 
by the production of a particle-antiparticle pair inside the potential barrier. It is different 
from the usual tunneling effect since it occurs inside the energy zone of the Klein paradox 
[59, 60]. Taking the quadratic form of the above equation for a generic scalar potential ^(a;),, 
we obtain the analogous Klein-Gordon equation, 

(ido - Vix)Y (f>{k, x) = {-d^ + m^) (f>{k, x), (89) 

which, from the mathematical point of view, due to the second-order spatial derivatives, 
has similar boundary conditions to those ones of the Schrodinger equation and leads to 
stationary wave solutions characterized by a relativistically modified dispersion relation. 

By depicting three potential regions by means of a rectangular potential barrier V{x), 
V{x) — Vo a < X < L, and V{x) = if x < and x > L, differently from the non- 
relativistic (Schrodinger) dynamics, we observe that the incident energy can be divided into 
three zones. The above-barrier energy zone, E > Vq + m, involves diffusion phenomena 
of oscillatory waves (particles). In the so called Klein zone [54, 59], E < V„ — m, we 
find oscillatory solutions (particles and antiparticles) in the barrier region. In this case, 
antiparticles see an opposite electrostatic potential to that seen by the particles and hence 
they will see a well potential where the particles see a barrier [61, 62]. The tunneling zone, 
Vo — m < E < Vq + m, ior which only evanescent waves exist [63, 64] in the barrier region, 
is that of interest in this work. In this context, our particular focus here is the calculation 
of phase times and dwell times and the possibility of a Hartman-like effect [10]. 

By evaluating the problem for this tunneling (evanescent) zone assuming that (f){k, x) are 
stationary wave solutions of the Eq. (89), when the peak of an incident (positive energy) 
wave packet reach the barrier x — at t — 0, we can usually write 

(l)i{k,x) = ex.p[i k x] + R{k, L) exp[—i k x] x < 0, 

(/>2(A;,x) = a{k) exp [—p{k) x] + /3{k) ex])\p{k) x] < x < L, 
(t)3{k,x) = T{k,L)ex.Y>[ik{x — L)] x > L, 

where the dispersion relations are modified with respect to the usual non-relativistic ones: 
k^^E^- and p{kf ^ - {E - K)'- 
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4>{k,x) 



The tunneling dynamics described by a relativistic wave equation can overcome all the 
(analytical) difficulties discussed in the last section. We demonstrate with complete math- 
ematical accuracy that, in some limiting cases of the relativistic (Klein-Gordon) tunneling 
phenomena where the relativistic kinetic energy is approximately equal to the potential en- 
ergy of the barrier, and mcL/h << 1, particles with mass m can pass through a potential 
barrier Vq of width L with transmission probability T approximately equal to one. Since 
T ~ 1, the analytical conditions for the stationary phase principle applicability which de- 
termines the tunneling (phase) time for the transmitted wave packets are totally recovered. 
Differently from the previous (non-relativistic) tunneling analysis, the original momentum is 
kept undistorted and there is no filter effect. We shall demonstrate that the tunneling time 
is then computed for a completely undistorted transmitted wave packet, which legitimizes 
any eventual accelerated transmission. 

A. Dynamics, variables and limits 

By assuming that the phase that characterizes the propagation varies smoothly around 
the maximum of g{k — ko), the stationary phase condition enables us to calculate the po- 
sition of the peak of the wave packet (highest probability region to find the propagating 
particle) . With regard to the standard one-way direction wave packet tunneling, for the set 
of stationary wave solutions given by Eq. (90), it is well-known [65] that the transmitted 
amplitude T{n,L) — \T{n,L)\ exp [i(/?(n, L)] is written in terms of 



for which we have made explicit the dependence on the barrier length L (parameter wL) 
and we have rewritten p{k) — wp{n), with p(n)^ = ^/T+'2r^^ — {rf — v/2). 

We illustrate the modified tunnehng transmission probabihties in the Fig. (8) for different 

propagation regimes {v = 0{NR), 1, 2, 5, 10)) by observing that the tunneling region is 
comprised by the interval {n — v/2y < 1, n > 0. 

The additional phase ip{n, L) that goes with the transmitted wave is utilized for calcu- 




(90) 



where we have suppressed from the notation the dependence on A;, and 




(91) 
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lating the transit time t^p of a transmitted wave packet when its peak emerges at x = L, 

^ dfe dn(fc) dy(n, £) ^ (L) 1 dy(n, £) 

~ d£;(fc) dA; dn ~ w{L) dn ' ^ 
evaluated at A; = /cq (the maximum of a generic symmetrical momentum distribution g{k—ko) 
that composes the incident wsive packet). By introducing the c/assica^ traversal time defined 

as r^^^ = L{dk/dE{k)) = L/v, we can obtain the normalized phase time, 

(93) 



_ f{n,L) 



where 

/(n, L) = 8n' [(2 + 8n'i; + v^) - (4n' + 3v) VTT^n^ 

+ 4 [(4 + in^v + v^) VlT^ -2v{2 + 3n^v)] Sh{p{n)wL) Ch{p{n)wL) 
'- V /J p{n)wL 

and 

g{n,L) = 16n' [2 (1 + 2n^v) - Vl + 2n^v (2n' + v)] 

+ 2 [(4 + 8n^v + v^) VI + 2n'v - 4x; (1 + 2n^v)\ Sh{p{n)wLf 

where Ch{x) = cosh {x) and Sh{x) = sinh (x). 

To illustrate the above results we plot the tunneling phase times in correspondence with 
the transmission probabihties of the Fig. (9) for the same different propagation regimes {v — 
0{NR), 1, 2, 5, 10)). 

By observing the results illustrated by the above figures, we can notice the possibility 
of accelerated (t^ < t^-^^-)), and eventually superluminal (negative tunneling delays, t^p < 0) 
transmissions without recurring to the usual analysis of the opaque limit {p{n)wL oo) 
which leads to the Hartman effect [10]. In the NR dynamics (Schrodinger equation solutions), 
the opaque limit and the consequent superluminal interpretation of the results of such an 
approximation (Hartman effect) were extensively discussed in the literature. Superluminal 
group velocities in connection with quantum (and classical) tunneling were predicted even 
on the basis of tunneling time definitions more general than the simple Wigner's phase time 
[19] (Olkhovsky et ai, for instance, discuss a simple way of understanding the problem 
[12]). Experiments with tunneling photons and evanescent electromagnetic waves [5, 22, 23] 
have generated a lot of discussions on relativistic causality, which, in addition to several 
analytical limitations, have ruined some possibilities of superluminal interpretation of the 
tunneling phenomena [13, 50, 65]. In a causa/ manner, the above arguments might consist 
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in explaining the superluminal phenomena during tunnehng as simply due to a reshaping of 
the pulse, with attenuation, as already attempted (at the classical limit) [49] , i. e. the later 
parts of an incoming pulse are preferentially attenuated, in such a way that the outcoming 
peak appears shifted towards earlier times even if it is nothing but a portion of the incident 
pulse forward tail [5, 50]. 

The Hartman effect is related to the fact that for opaque potential barriers the mean 
tunneling time does not depend on the barrier width, so that for large barriers the effective 
tunneling-velocity can become arbitrarily large, where it was found that the tunneling phase 
time was independent of the barrier width. It seems that the penetration time, needed to 
cross a portion of a barrier, in the case of a very long barrier starts to increase again after 
the plateau corresponding to infinite speed proportionally to the distance [82]. 

We do not intend to extend on the delicate question of whether superluminal group- 
velocities can sometimes imply superluminal signalling, a controversial subject which has 
been extensively explored in the hterature ([12] and references therein). Otherwise, the 
phase time calculation based on the relativistic dynamics introduced here offers distinct 
theoretical possibilities in a novel scenario, for the limit case where p{n) L tends to (with 
L ^ 0), in opposition to the opaque limit where p{n) tends to oo. Let us then separately 
expand the numerator f{n,L) and the denominator g{n,L) of the Eq. (93) in a power series 
of p(n) wL {p{n) — > 0) in order to observe that in the lower (upper) limit of the tunneling 
energy zone, where ri^ tends to i;/2 + (— )1, the numerical coefficient of the zero order term 
in p{n) wL amazingly vanish in the numerator as well as in the denominator! Since the 
coefficient of the linear term also is null, just the coefficient of the second order terms plays 
a relevant role in both series expansions. After expanding the Eq. (93), such a step-by-step 
mathematical exercise leads to 

4 [(4 + An'v + v') yi + 2n^v - 2v (2 + 3n'v)] , , , 

— = -7 : i 0(p(n)wL) (94) 

T^k) 3 [(4 + 8n^v + v^) VI + 2n^v - 4i; (1 + 2n^v)] ^^^^ ' ' ^ ' 

for small values of p{n). 

In parallel to the above way of calculating the traversal time, we believe that a more 

complete or, at least, complementary answer to such an inquiry is in the definition of the 

dwell time for the same tunneling configuration evaluated with a the Dirac equation in the 

Klein-Gordon equation form. The dwell time is a measure of the time spent by a particle in 
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the barrier region regardless of whether it is ultimately transmitted or reflected [26], 

tn = jjyx\Mk,x)\' (95) 

where jin is the flux of positive energy incident particles and 02 (fc, x) is the stationary state 
wave function inside the barrier. The normalized dwell time is thus given by 

to foijijL) 



where 



(96) 



p{n) J \ p{n) ) p{n)wL 
and 

Sh{pin)wLf 



9o{n,L) = 2Vl + 2n^ 



V 



1 + 



4:n'^p{nf 

We have computed the dwell times for the solutions of the relativistic wave equation (89) and 
we have illustrated it in the Fig. (10) in correspondence with the transmission probabilities 
of the Fig. (8). 

By taking the same limits that we have considered for the phase time, the expressions 
for the dwell time provide us with, 

lim — = l ^ / , ■ n'^t;/2^1, n\ v > 0. (97) 

from which we plot the comparative results in the Fig. (11). 

The transmission probability depends only weakly on the barrier height, approaching 
the perfect transparency for very high barriers, in stark contrast to the conventional, non- 
relativistic tunneling where T{n,L) exponentially decays with the increasing Vq. Such a 
relativistic effect is usually attributed that a sufficient strong potential which, being repulsive 
for electrons, is attractive for positrons, and results in positron states inside the barrier, 
which align the energy with the electron continuum outside [69] . 

B. Phase times and Dwell times - The relativistic case 

In order to reproduce the variational theorem for the Klein-Gordon equation (89), we 
could write 

{ido-E)(f>{k,x) = (98) 
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, its first derivative with respect to 



and, after some simple mathematical manipulations, its second derivative, 

+ 2E(f){k,x) = 0. 



dE 

By following the one- dimensional analysis here considered, it is easy to find that 

= 2{E - Fo)0V 



(99) 



(100) 



dE dx^ 



dx^ dE 



(101) 



dx \ dE dx ^ dEdx J 
where we clearly notice the presence of £^ — in place of m which appears in the non 
relativistic analysis [51] since in the strong non-relativistic limit E — Vq ^ m. Upon inte- 
gration over the length of the barrier we find 

dEdx) 



2{E-Vo) f dx\Mk, 
Jo 



x)\^ 



d(l)d(j^ ,+ d^ 

- r- 



(102) 



dE dx 

In front of the barrier [x < 0), the wave function consists of an incident and a refiected 

component given by 0i(fc,a;), and behind the barrier {x < L), there is only the transmitted 
wave 03(fc,.r) (see Eq. (90)). Under these conditions we evaluate the right-hand side of 
Eq. (102) as 



-2ik 



-^i\R\' + m + i[i\R\' + \T\ 



^dip{k,L) _^ Im[R] 



dk 



k 



dk_ 
dE' 



Since + \T\^ = 1, Eq. (102) becomes 



(E - Vo) £dx\Mk,x)\' = + ^Im[R] 



dip{k,L) ^ k 
' E' 



which gives 



dE 



t^f) _ tjf^ 1 Im[R] 
^■(fe) T-(fe) ''"(jk) E 



(103) 



(104) 



(105) 



where we have introduced the re-scaled dwell time, 

4"^ = ^^t^-) = rdxlUKxW (106) 

m k Jq 

which can be related to the correct definition of the probability density for the Klein-Gordon 
equation, 

jo = r dx [4{do<l>2) - {docf>l)<l>2]/ t^n^ = Uo/jin) = (jo/k), (107) 
Jo 
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and leads to the usual definition t^°^ in the non-relativistic hmit {^^^ i— > 1). In Eq. (106), 
the squared modulus of the wave function transforms as a Lorentz scalar, E—Vq transforms as 
a time-like component, and the integrand dx as well as k transform as space- like components 
of a Lorentz four-vector. It means that t^^"* has the correct Lorentz character since it 
transforms as a time-like component, which does not occur for t(^) of Eq. (95). 

As in the non-relativistic case [13, 51], the first term of Eq. (105) corresponds to the 
phase time or the aforementioned group delay. The second term comes from the explicit 
computation of the dwell time. However, the presence of the multiplicative factor 
in Eq. (106) introduces some novel aspects in the interpretation of the additional term 
—Im[R]/E as a self-interference term which comes from the momentary overlap of incident 
and reflected waves in front of the barrier [13]. As we can observe in the example illustrated 
in the Fig. 12, and by the usual definition (95), the dwell time is always positive. The 
re-scaled dwell time modulated by ( '^~^° changes sign when the total energy E equalizes the 
potential energy V^: an energy region comprised by the tunneling energy zone of the Klein- 
Gordon equation. Consequently, differently from the results we get from the non-relativistic 
analysis, the dwell time is not obtained from a simple subtraction of the supposed self- 
interference delay from the phase time that, in some circumstances [57, 65] describes 
the exact position of the peak of the scattered wave packets. 

At least for the moment, the above (relativistic) results do not necessarily demand for a 
confront with the (non-relativistic) predictions derived from the opaque limit analysis which 
results in the filter effect and the superluminal tunneling. To clear up this assertion, it 
is convenient to recover the limiting configurations (n^ v/2^1) of some of our previous 
results [57] for which the tunneling transmission probability (90) can be approximated by 



lim \T(n,L)\ 



1 + 



^ V » 1 

(wL) 



^ [l + (mLf]-K (108) 



from which, avoiding any kind of filter effect, we recover the probability of complete tunneling 

transmission when mL << 1, once we have \T{n,L)\ ^ 1. For the correspondent values of 

the phase times we obtain [57], 

t(.f) 4 1 

lim = - n'^v/2Tl, n\v>0, (109) 

that does not depend on mL, and we notice that its asymptotic (ultrarelativistic) limit always 
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converges to 0. Curiously, in the lower limit of the tunneling energy zone, — > t'/2 — 1, it 
is always negative. Since the result of Eq. (109) is exact, and we have accurately introduced 
the possibility of obtaining total transmission {transparent barrier), our result ratifies the 
possibility of accelerated transmission (positive time values), and consequently superluminal 
tunneling (negative time values), for relativistic particles when mL is sufficiently smaller than 
1 (=^ T fti 1). By observing that the barrier height has to be chosen such that one remains in 
the tunneling regime, it is notorious that the transmission probability depends only weakly 
on the barrier height, approaching the perfect transparency for very high barriers, in stark 
contrast to the conventional, non-relativistic tunneling where T{n,L) exponentially decays 
with the increasing Vq. Obviously, the above results correspond to a theoretical prediction, 
in certain sense, not so far from the experimental realization. The above condition should 
be naturally expected since we are simply assuming that the Compton wavelength {h/{rnc)) 
is much larger than the length L of the potential barrier that, in this case, becomes invisible 
for the tunneling particle. In general terms, the relativistic quantum mechanics establishes 
that if a wave packet is spread out over a distance d» 1/m, the contribution of momenta 
IpI ~ m >> l/d is heavily suppressed, and the negative energy components of the wave 
packet solution are negligible; the one-particle theory is then consistent. If we want to localize 
the wave packet in a region of space (wave packet width d) smaller than or of the same size 
as the Compton wavelenght, that is d < 1/m, the negative energy solutions (antiparticle 
states) start to play an appreciable role. The condition d < L < 1/m (where d < L is 
not mandatory) imposed over a positive energy component of the incident wave packet in 
the relativistic tunneling configuration excite the negative energy modes (antiparticlcs) and, 
qualitatively, report us to the Klein paradox and the creation of particle-antiparticles pairs 
during the scattering process which might create the intrinsic (polarization) mechanisms for 
accelerated and/or non-causal particle teletransportation. 



VII. SUMMARY AND OUTLOOK 



We have examined the meaning of the group delay in barrier tunneling for non-relativistic 

dynamics, extended to some relativistic dynamics. Our analysis was constructed upon the 
wave packet multiple peak decomposition which emerges when we study the one-dimensional 
potential scattering/ tunneling. Concentrated on some incongruities with the application of 
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the method of stationary phase for deriving phase times, we have extended the analysis to 
the study of the relativistic dynamics of the tunnehng phenomena. 

In the first part of our study, we have employed the SPM for analyzing the above-barrier 
diffusion problem. The main point we have elaborated on shows that the results of the 
SPM depend critically upon the manipulation of the amplitude prior to the application of 
the method. The method is inherently ambiguous unless we know, by some other means, 
at least the number of separate peaks involved. Essentially, what we have demonstrated is 
that the barrier results can be obtained by treating the barrier as a two-step process. This 
procedure involves multiple reflections at each step and predict the existence of multiple 
(infinite) outgoing peaks. 

In the second part of our analysis, wc have attempt to the physical framework where the 
study of a tunneling process is embedded. In the same analysis which generically involves 
analytically-continuous Gaussian pulses, holomorphic functions which do not have a well- 
defined front, we have quantified some aspects inherent to the quoted filter effect. In this 
case, we have observed that the amplitude of the transmitted wave is essentially composed 
by the plane wave components of the front tail of the incoming wave packet which reaches 
the first barrier interface before the peak arrival. As we have noticed, the filter effect by 
itself leads to ambiguities in the interpretation of the wave packet speed of propagation. As 
an alternative possibility to partially overcome these incompatibilities, we have claimed for 
the use of the multiple peak decomposition technique. Thus, in the same framework where 
we have treated the above-barrier diffusion problem, we have suggested a suitable way for 
comprehending the conservation of probabilities for a very particular tunneling configuration 
where the asymmetry presented in the previous case was eliminated, and the phase times 
could be accurately calculated. An example for which, we believe, we have provided a simple 
but convincing resolution which may either reignite the tunneling time controversy or help 
settle the issue. 

In more general lines, there have also been some trying of yielding complex time delays 
for tunneling analysis, ultimately due to a complex propagation constant. In this framework, 
the supposition of superluminal features is considered artificial since the transmitted peak is 
not related causally to the corresponding incoming peak [3]. In parallel to the most sensible 
candidate for tunneling times [1, 3], a phase-space approach have been use to determine 
a semi-classical traversal time [66-68]. This semi-classical method makes use of complex 
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trajectories which, by its turn, enables the definition of real traversal times in the complex- 
ified phase space. It is also commonly quoted in the context of testing different theories for 
temporal quantities such as arrival, dwell and delay times [1, 3] and the asymptotic behavior 
at long times [11, 14]. 

As a not so ambitious but suggestive possibility to partially overcome some of the incon- 
gruities here noticed and quantified, which appear when we discuss tunneling phase times, 
we have claimed for the use of the multiple peak decomposition [16] technique developed for 
the above-barrier diffusion problem. Thus we have introduced a suitable way for compre- 
hending the conservation of probabilities for a very particular tunneling configuration where 
the asymmetry presented in the previous case was eliminated, and the phase times could be 
accurately calculated. We let for a subsequent analysis the suggestive possibility of inves- 
tigating the validity of our approach when confronting with the intriguing case of multiple 
opaque barriers [46], in particular, in the case of non-resonant tunneling. Still concerning 
with the future theoretical perspectives, the symmetrical colliding configuration also offers 
the possibility of exploring some applications involving soliton structures. In summary, the 
above discussion reinforces the assertion that the investigation of wave propagation across 
a tunnel barrier has always been an intriguing subject which is wide open both from a 
theoretical and an experimental point of view. 

We have concluded that the transmission probability depends only weakly on the barrier 
height, approaching the perfect transparency for very high barriers, in stark contrast to the 
conventional, non-relativistic tunneling where T{n, L) exponentially decays with the increas- 
ing Vq. Such a relativistic effect is usually attributed to a sufficient strong potential that, 
being repulsive for electrons, is attractive for positrons and results in positron states inside 
the barrier, which align the energy with the electron continuum outside [69]. 

Finally, in the last part of the manuscript, we have been concentrated on the accurate 
calculation of phase times and dwell times for the energy zone of tunneling governed by a rel- 
ativistic (Dirac/Klein-Gordon) wave equation. The well-known Klein paradox - unimpeded 
penetration of relativistic particles through high and wide potential barriers - is one of the 
most exotic and counter intuitive consequences of quantum electrodynamics (QED). The 
phenomenon is discussed in many contexts in particle, nuclear and astrophysics but direct 
tests of the Klein paradox using elementary particles have so far proved impossible. Here we 
have quantified the conditions for the occurrence of accelerated and, eventually, superluminal 
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tunneling transmission probabilities in order to eliminate the called filter (Hartman) effect 
from the transmitted waves. By eliminating the filter effect, the transmission probabilities 
approximates the unitary modulus (complete transmission through a transparent medium) . 
As a result, the calculation of the phase time could be accurately evaluated in order to give 
the exact position of the transmitted wave packet (or particles), for which, in the standard 
non-relativistic analysis, due to the distortion of the original momentum distribution, the 
position of the peak is shifted. 

By considering the magnitude of the parameter mL {mc^ /[h{c/ L)] in standard units) 
for an electron with mass ~ 0.5 MeT^, and observing that in natural units we have 
Q.2MeV pm ~ 1, we conclude that it should be necessary a potential barrier of width 
L « 1pm to permit the observation of the quoted superluminal transmission. By princi- 
ple, its observation makes the effect relevant only for some exotic situations as, for instance, 
positron production around super-heavy nuclei {Z ~ 170) [70] or evaporation of black holes 
through generation of particle- antiparticle pairs near the event horizon [71]. The super- 
position of several barriers in a nanoscopic scale could perhaps provide a truly observable 
effect. 

In the most common sense, the above condition should be naturally expected since we 
are simply assuming that the Compton wavelength {h/ (mc)) is much larger than the length 
L of the potential barrier that, in this limit situation, becomes invisible for the tunneling 
particle. The relativistic quantum mechanics establishes that if a wave packet is spread 
out over a distance d » 1/m, the contribution of momenta \p\ ~ m » l/dis heavily 
suppressed, and the negative energy components of the wave packet solution are negligible; 
the one-particle theory is then consistent. However, if we want to localize the wave packet 
in a region of space (wave packet width d) smaller than or of the same size as the Comp- 
ton wavelenght, that is o? < 1/m, the negative energy solutions (positron states) start to 
play an appreciable role. This qualitative arguments report us to the Klein paradox and 
the creation of particle- antiparticles pairs during the scattering process which might create 
the intrinsic (polarization) mechanisms for accelerated and/or non-causal fermion teletrans- 
portation. The condition d < L < 1/m imposed over a positive energy component of the 
incident wave packet in the relativistic tunneling configuration excite the negative energy 
modes (antiparticles) in the same way that the motion of electrons in a semi-conductor is 
concatenated with the motion of positively charged holes. 
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Turning back to the context of the nanoscopic scale structures, the most challenging 
possibility of observing similar effects occurs for massless (or effective mass) Dirac fermions 
in graphene structures. Even though the linear spectrum of fermions in graphene implies 
zero rest mass, their cyclotron mass approaches to lO^^me [72], which increases the super- 
luminal tunneling scale to 1 angstrom. In spite of the theoretical focus, the results here 
obtained apply to some configurations which should deserve further attention by exper- 
imenters in the study of the graphene structures where the dynamics of the electron is 
described by a relativistic-like dynamics. For bi-layer structures, due to the chiral nature 
of their quasiparticles, quantum tunneling in these materials becomes highly anisotropic, 
qualitatively different from the case of normal, non-relativistic electrons. In the last years, 
it has been speculated that, from the experimental point of view, the graphene provides 
an effective medium for mimicking relativistic quantum effects where, for instance, mass- 
less Dirac fermions allow a close realization of Kleins gedanken experiment whereas massive 
chiral fermions in bilayer graphene offer an interesting complementary system that eluci- 
dates the basic physics involved. In conventional two-dimensional systems, strong enough 
disorder results in electronic states that are separated by barriers with exponentially small 
transmittance [69]. In contrast, in single- and bi-layer graphene materials all potential bar- 
riers are relatively transparent {T{n,L) fa 1): the quasiparticles in graphene exhibit a linear 
dispersion relation E — hkvf that corresponds to the pseudo-ultrarelativistic limit of our 
analysis for pseudo- massless particles traveling with their Fermi velocity Vf. In this case 
there are pronounced transmission resonances where T approaches unity for some particular 
geometric configurations, which does not allow charge carriers to be confined by potential 
barriers that are smooth on atomic scale. Some authors have demonstrated experimentally 
and theoretically, that the biased graphene bilayer is a tunable semiconductor where the 
electronic gap can be controlled by the electric field effect reaching values as large as 0.3 eV 
[73]. These results can have important implications for the development of carbon-based 
electronics. 

Once it has been observed that the physical essence of the Klein paradox hes in the predic- 
tion that particles can pass through large repulsive potentials without exponential damping 
[54, 55], in this scenario, the Klein-paradox and the accelerated tunneling transmission as- 
sociated with relativistic-like phenomena at nanoscopic scale can be tested experimentally 
using graphene devices. To sum up, as future perspectives we intend to investigate the ap- 
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pearance of an equivalent smaller effective mass value Mg// « m due the minimal coupling 
of the charged particle magnetic momentum with an external magnetic field, in particular, 
for the massless electron propagation in a single-layer graphene which could introduce some 
novel ingredients for quantifying these peculiarities of the relativistic tunneling effect. 
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FIG. 1: The phase derivative 8'(fc) dependence on k/w. Q'{k) does not present an adequate analytical 
behavior (smoothness) for the applicability of the SPM when k approximates to w since the phase derivative 
oscillates too rapidly (The phase is not stationary) . The method can be accurately applied for larger values 
of k/w when the phase is really stationary. 



50 



1. 

0.8 
0.6 
0.4 
0.2 






1. 




0.8 




0.6 








0.2 




Q 




1. 




0.8 




0.6 




0.4 


>> 


0.2 







?! 








1 


1. 






0.8 




0.6 






0.4 




0.2 








1. 

0.8 
0.6 
0.4 
0.2 

1. 
0.8 
0.6 
0.4 
0.2 

1. 
0.8 
0.6 
0.4 
0.2 



Mod(— —Y\ 




Mod( ■ f 

t(xo,0) 



Mod( f 

t(x„,0) 



Mod( 



"fix,!;) 

■9(xo,0) 



Modi ■ ? 

1l(x„,0) 



Uod( f 

1l(.x„,0) 



'ii.xo, 0) 









A 







-10 



10 



15 



0.8 
0.6 
0.4 
0.2 


1. 

0.8 
0.6 
0.4 
0.2 


1. 

0.8 

0.6 

0.4 

0.2 



1. 

0.8 

0.6 

0.4 

0.2 



1. 

0.8 

0.6 

0.4 

0.2 

1. 

0.8 

0.6 

0.4 

0.2 

1. 

0.8 

0.6 

0.4 

0.2 





x/a 



FIG. 2: Erroneous interpretation of the scattering of an incoming wave packet by a one-dimensional potential 
barrier of height Vq and width L. It has the purpose of illustrating the deficiencies inherent to the wrong 
applicability of the SPM. We have plotted the propagating wave packets at the corresponding times tn = 
{ma'^)[n {L/a)]/{aqo) (with n = 0, 1, 5 and with the normalization constraint mo^ = 1) by assuming the 
incoming wave packet starts at a; = —{ko L)/ (2go)- From the false behavior of the density of probabilities it 
becomes obvious that the total probability is not conserved as it was expected. The square of the amplitude 
modulus would supposedly represent a collision of a wave packet of average width a with a potential barrier 
Vo of width L = 5a where, for illustrating reasons, we have adopted ko = \/2w and wa = 10000. Only a 
fixed region in x close to the barrier is shown. 
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FIG. 3: Multiple peak decomposition of a propagating wave packet in the above potential barrier scattering 
problem. We have plotted the first few reflected and transmitted waves for times f„ = {ma'^){n{L/a)\/{aqo) 
(with n = 0, 1, 5 and ma^ = 1) in correspondence with Eq. (41), by assuming that the incoming wave 
packet starts at x = — (fco L)/{2qo). The density of probabilities represents the collision of a wave packet of 
average width a with a potential barrier of width L = 5a. Just for illustrating reasons, we have adopted 
ko = (■v/T0w)/3 {wa = 10000) and we have printed the wave packet amplification multiplying factor (A) 
(individually adopted for visual convenience for each wave packet) when necessary. 
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FIG. 4: The illustrative confront between the analytical (line) and the numerical (star symbols) results for 
reflected and transmitted waves when the condition 1 < ka/qo < -k is satisfled and, consequently, both the 
analytical and the numerical results have to be equivalent. The stars are from the numerical convolution 
of the plane- wave solution and the curves are from the analytical expressions for the first three i?„ and r„. 
The density of probabilities represents the collision of a wave packet of average width a with a potential 
barrier Vq of width L = 0.8a. Again for illustrating reasons, we have adopted kg = (5\/2w)/7 and we have 
printed the wave packet amplification multiplying factor (^4) when necessary. The reason for choosing a so 
large value for the parameter wa {wa = 10000) in all these figures is to guarantee the convergence between 
the analytical and the numerical results for all the decomposed peaks. 
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FIG. 5: Dependence of the wave packet shape on the cut ojf value of a momentum distribution 
centered around ko = 0.5w with the values of k comprised between and k„,t ^g. 
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FIG. 6: Normahzed phase times for a symmetrical ti?linehng configuration of the symmetrized wave 
function representing the colhsion of two identical bosons (black lines) and the antisymmetrized 
wave function representing the collision of two identical fermions (red lines) . The results are plotted 
in comparison with the standard one-way direction tunneling phase times (blue lines). These times 
can be understood as transit times in the units of the classical traversal time = {mL)/k. All 
the above phase time definitions present the same asymptotic behavior. 
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FIG. 7: Exact phase time (solid line), self-interfeftnce delay (dotted line), and the dwell time 
(dashed line) as a function of the normalized energy n = jvip' oc -Eo/K for the two identical 
particles (black line) and the standard one-way direction (red line) collision with a rectangular 
potential barrier. For the colliding two identical fermions (first plot) we have assumed that the 
wave function is totally antisymmetrized. For the colliding two identical bosons (second plot) we 
have assumed that the wave function is totally symmetrized. These times are normalized by the 
classical traversal time = {mL)/k, and here we have adopted wL = Att for a = wL\/\ — n. 
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FIG. 8: Tunneling TRANSMISSION probabilities |r(n,L)|^ for the incidence on a retangular po- 
tential barrier of height Vq for the dynamics of the relativistic wave equation versus the normalized 
energy n^{k). We have classified the energy zones by the line thickness: the thick line corresponds 
to the tunneling zone, the intermediate line corresponds to the Klein zone and the thin line cor- 
responds to the above-barrier zone. Here we have adopted the illustratively convenient value of 
wL = 2iT and have set u = 0, 1, 2, 5, 10, where NR regime can be parameterized by f = 0. We 
have constrained our analysis to n^(A:) > since we have assumed Vq > 0. 
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Normalized Energy n\k) = k^/ 2rnV 

FIG. 9: The related tunneling PHASE times in correspondence with their the tunneling TRANS- 
MISSION probability for the dynamics of the Klein-Gordon equation form. Again by means of 
the line thickness we observe that the tunneling region is comprised by the interval (n^ — < 
1, n > 0. We have used the same criteria and the same set of parameters from the Fig. (8). 
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FIG. 10: The related tunneling PHASE times in correspondence with their tunneling TRANSMIS- 
SION probability for the dynamics of the Klein-Gordon equation form. Again, by means of the line 
thickness we observe that the tunneling region is comprised by the interval (n^ — 1;/2)^ < 1, > 0. 
We have used the same criteria and the same set of parameters from the Fig. (8). 
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Phase Time n = u/2 +1 
Dwell Timen = u/2 +1 
Pfme Time n = u/2 - 1 
Dwell Time n = v/ 2 -1 




Nonmlized Energy n (k) = J^/lmV 



FIG. 11: Phase times (solid lines) and dwell times (dashed lines) in the limits of the tunneling 
energy zone as a function of the propagation regime (normalized energy parameter n^(A;)). The 
analytical transition between the tunneling zone and the above-barrier zone is given by = v/2+1, 
and between the tunneling zone and the Klein zone by = i;/2 — 1. In the lower limit of the 
tunneling energy zone, ra^ ^ f /2 — 1, it is always negative. 
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FIG. 12: Delay times calculated from the dynamics of the Klein-Gordon equation: Phase Time 
(Dash-dotted line), Dwell Time (Dashed black line), Self-Interference term (Solid line), and the 
re-scaled Dwell Time (Dashed red line). In fact the tunneling region is comprised by the interval 
(n^ — v/2y^ < 1, > 0. Here we have adopted the illustratively convenient value of wL = 27r with 
V = 5, in comparison with the non-relativistic results parameterized by u — > 0. . 
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TABLE I: The values of k numerically obtained by increasing the barrier extension 
L. The values are calculated in terms of the wave packet width a for different values 
of the potential barrier high expressed in terms of w a. We have fixed the incoming 
momentum by setting ko a equal to 1. 



w a 
L/a 




1.5 




2.0 




4.0 




6.0 




8.0 




10 




20 


0.00 


1 


0000 


1 


0000 


1 


0000 


1 


0000 


1 


0000 


1 


0000 


1 


0000 


0.05 


1 


0062 


1 


0188 


1 


1777 


1 


4156 


1 


6238 


1 


7726 


1 


9834 


0.10 


1 


0235 


1 


0648 


1 


3799 


1 


6769 


1 


8547 


1 


9397 


2 


0051 


0.15 


1 


0489 


1 


1223 


1 


5349 


1 


8251 


1 


9505 


1 


9937 


2 


0133 


0.20 


1 


0794 


1 


1825 


1 


6571 


1 


9178 


2 


0000 


2 


0204 


2 


0203 


0.25 


1 


1129 


1 


2420 


1 


7575 


1 


9813 


2 


0317 


2 


0390 


2 


0272 


0.30 


1 


1478 


1 


3001 


1 


8430 


2 


0289 


2 


0562 


2 


0551 


2 


0342 


0.35 


1 


1836 


1 


3565 


1 


9185 


2 


0679 


2 


0779 


2 


0704 


2 


0413 


0.40 


1 


2196 


1 


4116 


1 


9874 


2 


1025 


2 


0986 


2 


0857 


2 


0484 


0.45 


1 


2558 


1 


4657 


2 


0524 


2 


1350 


2 


1191 


2 


1012 


2 


0556 


0.50 


1 


2921 


1 


5194 


2 


1155 


2 


1668 


2 


1399 


2 


1170 


2 


0628 


0.55 


1 


3285 


1 


5729 


2 


1785 


2 


1988 


2 


1611 


2 


1331 


2 


0701 


0.60 


1 


3649 


1 


6266 


2 


2429 


2 


2314 


2 


1828 


2 


1495 


2 


0775 


0.65 


1 


4015 


1 


6809 


2 


3101 


2 


2651 


2 


2051 


2 


1663 


2 


0850 


0.70 


1 


4383 


1 


7360 


2 


3819 


2 


3002 


2 


2281 


2 


1834 


2 


0925 


0.75 


1 


4751 


1 


7920 


2 


4599 


2 


3367 


2 


2518 


2 


2009 


2 


1001 


0.80 






1 


8489 


2 


5466 


2 


3751 


2 


2761 


2 


2188 


2 


1078 


0.85 




* 


1 


9065 


2 


6456 


2 


4154 


2 


3013 


2 


2371 


2 


1155 


0.90 




* 


1 


9646 


2 


7627 


2 


4578 


2 


3272 


2 


2558 


2 


1234 


0.95 




* 




* 


2 


9091 


2 


5028 


2 


3540 


2 


2750 


2 


1313 


1.00 




* 






3 


1137 


2 


5501 


2 


3818 


2 


2917 


2 


1392 



Tor the values of L marked with *, we can demonstrate by means of Eqs. (61-62) that the 
modulated momentum distribution has already been completely distorted and the maximum 
loses its meaning in the context of applicability of the method of stationary phase. 
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